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A quenched second order phase transition is modeled by an effective <E> 4 -theory with a time- 
dependent Hamiltonian H(t), whose symmetry is broken spontaneously in time. The quantum field 
evolves out of equilibrium (nonequilibrium) during the phase transition as the density operator 
significantly deviates from p(t) = e _/3H ' i) /Zh- The recently developed Liouville-von Neumann 
(LvN) method provides various quantum states for the phase transition in terms of a complex 
solution to the mean-field equation, which is equivalent to the Gaussian effective potential in the 
< ^^> static case and to the time-dependent Hartree-Fock equation in the nonequilibrium case. Using 

the multiple scale perturbation theory (MSPT) we solve analytically the mean-field equation to the 
(3JT)' first order of coupling constant and find the quantum states during the quenched second order phase 

^ , transition. We propose a renormalization scheme during the process of phase transition to regularize 

the divergences, which originate from the mode coupling between hard and hard modes or between 
the soft and hard modes. The effect of mode coupling is discussed. 
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\ The dynamics of nonequilibrium (out of equilibrium) quantum fields has been a subject of much discussion in the 
recent years [Q,^). Nonequilibrium quantum fields play important roles in a variety of different scenarios. During a 
second order phase transition, the time scale of relaxation of the scalar field lags behind the time scale of the change of 
the effective potential. Consequently the field evolves out of equilibrium as it tries to relax to the new vacuum state, 
having a nonvanishing vacuum expectation value. Such nonequilibrium effects play a crucial role in topological defect 
Oh, formation both in condensed matter systems as well as in the early universe In a seminal paper, Kibble first 

showed how the correlation length is crucial in determining the initial density of topological defects j(| . These ideas 
were applied by Zurek who proposed that it may be possible to quantitatively test the Kibble mechanism of defect 
,£2 • formation in condensed matter systems like superfluid He 4 0. He argued that, due to the phenomenon of critical 
slowing down near the phase transition point, the correlation length relevant for determining the initial density of 
the defects is not the equilibrium correlation length at the Ginzburg temperature but that at the time when the field 
dynamics essentially freezes out. He also found a power law behaviour in the dependence of the correlation length 
(and consequently initial defect density) on the quench rate. There have been various attempts to experimentally test 
the Kibble-Zurek prediction of initial defect density ||. A comprehensive review of these issues is given in Ref. 

The nonequilibrium evolution of a scalar field also is responsible for the establishment of large-scale correlations 
leading to growth of domains ||. This is particularly relevant in the context of Disoriented Chiral Condensate (DCC) 
formation, where only a large domain of DCC can produce the dramatic fluctuations in the charged to neutral pion ratio 
p0[ . The field theoretical treatment from nonequilibrium to equilibrium is particularly important in understanding 
how a system of quarks and gluons thermalizes after a heavy-ion collision to form a new state of matter called quark- 
gluon plasma. Searches of quark-gluon plasma are currently underway at RHIC in Brookhaven and LHC in CERN. 
In another context, the recent preheating mechanism of the Universe after inflation might have accompanied the 
explosive decay of the infiaton into a large number of soft quanta of the field(s) [O-O. This process, leading to 
a large population of the low momentum modes, occurs out of equilibrium. A subsequent interaction between the 
various modes leads to the transfer of energy from low to high momentum modes and the eventual thermalization of 
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the Universe. Another exciting aspect of the nonequilibrium field dynamics which has recently received much attention 
is the issue of decoherence and consequent appearance of classicality during a nonequilibrium phase transition [fl5|-|l7| . 

A feature common to all such (quenched) phase transitions is the evolution far away from equilibrium and the 
exponential growing of the soft (long wavelength) modes, which necessarily leads to the domain growth. Finite 
temperature field theory based on equilibrium or quasi-equilibrium methods does not describe all the processes of 
nonequilibrium evolution even though the imaginary part of the complex effective action yields the decay rate [ ^8| . 
To treat properly such nonequilibrium quantum evolution, Schwinger and Keldysh first introduced the closed-time 
formalism p9| . The closed-time formalism or the 1/N expansion method have been applied to the nonequilibrium 
<fr 4 -theory to explain the phenomenon of domain growth ||^0) . In a recent paper Q] , a canonical method called the 
Liouville-von Neumann (LvN) approach has been developed that unifies both the functional Schrodinger equation for 
quantum evolution and the quantum LvN equation for quantum statistics. At the lowest order of coupling parameter 
the LvN approach yields the Gaussian approximation for static quantum fields and leads to mean-field equations for 
nonequilibrium quantum fields that can be cquivalcntly obtained by the time-dependent Hartree-Fock method. 

The purpose of this paper to study the renormalization of an effective theory of <I> 4 field undergoing a quenched 
second order phase transition. The $ may represent a quantum field of an effective system interacting with an 
environment or an order parameter of the system. In all such cases our model for the second phase transition is 
described by the Hamiltonian 



The transition is modeled through a sudden and instantaneous change in sign of m| at t — 0; this change destabilizes 
the field whose true vacuum was at $ = before t < and signifies the onset of the spontaneous symmetry breaking 
transition. Before the phase transition the quantum field is in a static state of either vacuum state or thermal 
equilibrium and can be accurately described by the Gaussian effective potential (GEP) |2^,^]. But the onset of 
the quench drives the field out of equilibrium, so the finite temperature field theory or GEP may not be directly 
applied; instead the LvN or the Hartree-Fock method provides the correct evolution of the field even after the phase 
transition. In this paper we shall employ the LvN method, as this method exactly recovers the GEP result before the 
phase transition and the time-dependent Hartree-Fock result after the phase transition. 

The Hartree-Fock method has been popular and useful in studying the nonequilibrium evolution of the field |]9|, po| , p^| . 
Even though the renormalization is well understood even for nonequilibrium fields in Refs. P j2(i| , p^ |, a systematic and 
explicit renormalization scheme has not been properly addressed for the phase transitions. In this paper we propose 
such a renormalization scheme for phase transitions by applying the multiple scale perturbation theory (MSPT) 
method |^,^6| to the LvN method and thus obtain the renormalized parameters after phase transitions. The MSPT 
method introduces the multiple time scales determined by the nonlinear coupling constant of the system, say A, 
and leads to the solution with the correct renormalized frequency even at order Xh. In field theoretical terms the 
MSPT solution even at the lowest order is equivalent to summing selected diagrams to all orders and is an accurate 
approximation for a sufficiently long time, l/(Xh). The MSPT method turns out particularly useful within the LvN 
method since the procedure of renormalization in the mean-field or the Hartree-Fock equation amounts to obtaining 
the renormalized frequencies of the hard and soft modes according to the MSPT method to zeroth order. The 
calculation is carried out within the Gaussian approximation and, therefore, does not address all the issues pertaining 
to the non-Gaussian effect of nonlinear mode mixing. It nevertheless provides an analytical (albeit, perturbative) 
method of obtaining expressions for the renormalized parameters during a spontaneous symmetry breaking phase 
transition. 

The paper is organized as follows. In Sec. II, we show how the LvN method can be used to obtain the renormalized 
Gaussian effective potential. We also obtain the appropriate form of the renormalized mass and coupling constant 
using the GEP. In Sec. Ill, we study a simple quantum-mechanical model of two coupled quartic oscillators as a 
precursor to studying the nonequilibrium evolution of the quantum fields during a quenched second order phase 
transition. We apply the MSPT method to obtain the correct frequency for each mode. Although the issue of 
renormalization is not relevant for this quantum mechanical model, it nevertheless provides some insight into the 
nonlinear effect of coupling between soft and hard modes during the quenched second order phase transition. In Sec. 
IV, we deal with the renormalization of mass and coupling parameters during the nonequilibrium evolution of the 
quantum fields without phase transitions. In Sec. V, by applying the MSPT we obtain the renormalized solution to 
field equations after the second order phase transition. A summary and discussion of our results is given in Sec. V. 




(1) 



where 




t < 0, 
t > 0. 
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II. GEP FROM THE LVN METHOD 



The essential idea of the LvN method plj is to simultaneously solve the functional Schrodinger equation 

d 



ih m \*(t)) = H(t)m)h 



(2) 



and the quantum LvN equation 



ih-0(t) + [0(t),H{t)] 



(3) 



The LvN method is based on the observation by Lewis and Riesenfeld [^?J that the eigenstates of 0(t) are the 
exact quantum states of time-dependent Schrodinger equation (^) up to time-dependent phase factors. For a time- 
independent system the Hamiltonian itself satisfies the LvN equation and its energy eigenstates are exact quantum 
states as expected. In terms of appropriately selected operators satisfying Eq. one may find the density operator 
and the Hilbert space of exact quantum states. In this sense the LvN method unifies quantum statistical mechanics 
with quantum mechanics. The LvN method treats the time-dependent, nonequilibrium, system exactly in the same 
way as the time-independent, equilibrium, one. Also the LvN method can be applied to nonequilibrium fermion 
systems with a minimal modification p8| . 

To compare the LvN method with the well-known GEP [p3| , a nonperturbative method, we consider a field theoretic 
model to which both methods can apply. Such a model is provided by the <I> 4 -theory with the Hamiltonian 



H(t) 



d 3 x 



(2tt 



,3/2 



^$4 

4! 



(4) 



In the GEP method the effective potential of a classical background field <j) c includes corrections from quantum 
fluctuations. The field is assumed to fluctuate in a Gaussian vacuum state with the expectation value 

(o|$|o)q = ^ c , <o|n|o) G = o, (5) 

where the field is divided into the classical background field and fluctuations 

*(t,x) = & + x). (6) 
To express Eq. ^ as an infinite sum of coupled quartic oscillators, we redefine the Fourier modes of fluctuation as 

4 +) (*) = + 0-k(t)], = 0-k(*)], (7) 

where 

d 3 k 



(2tt)3/ : 



Then the mode-decomposed Hamiltonian is given by 

1 o 



a 

where a = {(±),k} and 



Ab 
4! 



B i2 



4f^ 



2 



(8) 



(9) 



(10) 



Here the odd power terms of <j) a are neglected since we are interested in symmetric quantum states. Following Kim 
and Lee we introduce the operators 



Al{t) = --j=[(p a (t)7t a - <p a (t)$ a ], A a (t) = -l=[^* a (t)fr a - <j>Z(t)4> a ], 



(11) 
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with overdots denoting time derivative, and express the Hamiltonian (Q) as normal-ordered quadratic and quartic 
parts in A and A': 



where 



H = Hq -\- Hp 



Hg = g + ul^)Al + + ^M(2AiA a + 1) + + 

a. 

EtofrMM + ^M2AiA a + 1) + < 2 it 2 ] + nA € + Ab ^ 



A R ft 2 



a,/3 



(12) 



(13) 



and 



, Q fc=Q ^ ' 



4! 



(14) 



The strategy of the GEP and LvN methods is to exactly solve the quadratic part Hq and perturbatively treat the 
quartic part Hp. In the GEP method, the GEP is obtained by minimizing the free energy with respect to a variable 
frequency, which consequently determines ip a . In the LvN method, the requirement that the creation and annihilation 
operators ( pi] ) satisfy the LvN equation (|^) leads to the mean-field equations for ip a . 

In the GEP method for a static system with time-independent coupling parameters ]22] , p3| ], one minimizes the 
vacuum energy with respect to a trial Gaussian state of the form 



<>a) 



/ \ 1/4 




V n J 




~2h^ 



(15) 



where Q a are free parameters. The variational parameters fl a are time-independent for the static system with time- 
independent coupling parameters, but depend on time for a noncquilibrium system with time-dependent coupling 
parameters. The time-dependent variational principle by Dirac P9[ is required for the time-dependent system, to 
which we shall apply the LvN method in this paper. The time-dependent phase factor in Eq. (US) has been inserted 



to satisfy the time-dependent Schrodinger equation. Indeed, the Gaussian state 
by A in Eq. dill), which can be shown by choosing the solution 



1 



-iQ a t 



is the vacuum state annihilated 



(16) 



and writing Eq. (|iT|) in the canonical form 



4w 



-iQ a t 



A a (t) 



AQ. a t 



2hfi a 



(17) 



Then the Hamiltonian leads to the Gaussian effective potential 



Ab^ 2 



B aA 



4^1 ' Q 8 \^ 2il a 2 

a \ / \ a / 

The minimization of the Gaussian effective potential with respect to £l a leads to the gap equation 



(18) 
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where 



Ate) 



^ B J.2 



'B 



A B fi 
2 



(19) 



(20) 



Now the Gaussian effective potential takes the form 

.,2 



A B ?i 2 



where 



'in 



l-2n 



d 3 k 



(2tt) 3 / 2 O^- 2 ™ ' 



(21) 



(22) 



In the GEP method the renormalization of coupling parameters (constants) are prescribed by matching the coeffi- 
cients of the effective potential with renormalized ones. Following Ref. ^3|, the renormalized mass in the Gaussian 
state, a symmetric state around <f> c = 0, is given by 



d 2 V G 



■ a =o 

2 , AbTi t , , 

Similarly, the renormalized coupling constant is found to be given by 

d 4 V G 



(23) 



Ar = 



[1 + Y J -iO*)]" 



(24) 



The case of particularly relevance to our model with a spontaneous symmetry breaking is when Ab is positive and 
finite [^3| . The other cases corresponding to different values of Ab have also been discussed in detail pj| . In our model 
the time-dependent mass term spontaneously breaks the symmetry during the quenched phase transition, which may 
be regarded as a consequence of the interaction with an environment. 

On the other hand, in the LvN method, the operators ( |li"|) are required to satisfy the LvN equation (J3[) for the 
truncated Hamiltonian Hq. This leads to the mean- field equation for each mode 



Ab^- 







If a = 0. 



(25) 



The LvN method allows even time-dependent coupling parameters as long as Eq. ( |25|) is satisfied. Further, the 
operators (pd|) are endowed with the standard commutators 

(26) 



[l a (i),4(i)] =8 aJl , 
which are guaranteed by the Wronskian conditions 



(27) 



The Fock space for each mode is constructed from the time-dependent creation and annihilation operators, A^it) and 
A a (t). The Gaussian vacuum state of ath mode, which is an approximation to the true vacuum state, is annihilated 
by !„(*): 
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A a (t)\O a ,t) G = 0, 



(28) 



and has the coordinate representation 



U/4 











— ex P 


2h ^ 



Similarly, the nth excited state is obtained by applying A^ a (t) n times: 

\n a , t) = n L= Ui (t)) na \o ai t) 

and, in the coordinate representation, is given by 

1 1 



n+l/2 



l>a,t) 



((2ft)«n!)V 2 {2*1uplip a yt*W a 







j exp 


2h V * a 9a _ 



(29) 



(30) 



(31) 



where H n is the Hermite polynomial. 

A few comments are in order. The quantum states ( ^9|) and ( |3l|) in general depend on time, except for trivial 
time-dependent phase factors, whenever the coupling parameter tob or Ab depends on time. In fact, these states 
constitute the time-dependent Fock space. The Gaussian vacuum state of the field <E> is the product of Gaussian state 
of each mode: 



io,t) G -ni°«' t ) ( 



(32) 



Similarly, the excited state of the field is the product of each mode state, at least one state being excited. This Fock 
representation of quantum field is unitarily inequivalent because all excited states of the field are not countable |30) . 
Further, there is another feature of the noncquilibrium field that is absent from the static quantum field. In spite 
of the fact that each Fock space at two different times is unitary equivalent through a Bogoliubov transformation, 
the Gaussian vacuum states of the field are orthogonal to each other at two different times due to infinite number of 
modes: 



! (0,i'|0,t) G =0 ) (t'^t). 



(33) 



Now the task in the LvN method is to solve the mean-field equation (g5|). For the static system the complex function 
( |l6| ) is a solution to Eq. (|2^) if fl a satisfies the gap equation ( |l9] ) . This shows the equivalence between the LvN and 
GEP methods for at least static case. Also it follows that the mean-field equations ( p5| ) for symmetric quantum states 
with cj) c = 0, upon mass renormalization (J23|), become renormalized 



where the renormalized frequency is 



o 2 — o 2 



Ra 



0. 



(34) 



(35) 



Therefore, a simple renormalized solution is given by 



V2TT 



Rq 



(36) 



This implies that the renormalized Gaussian vacuum state (^9|) and its excited states (^Tj) for the static Hamiltonian 
(|J) indeed describe the free theory with the renormalized mass. 

We now compare the Gaussian effective potential ( p^ ) with the conventional perturbation theory. This comparison 
provides a direct interpretation of the results of the MSPT method, which coincide with the results of GEP method 
to first order of Xb^, as will be shown in Sees. IV and V. The renormalized mass is particularly of interest from the 
point of view of instability during the phase transition. The conventional perturbation theory is based on the solution 



G 



It is a solution to Eq. ( |25| ) without the self-interaction term (Ab = 0) and leads to a different Gaussian vacuum state 
(p9|), which is denoted as |0, i)my So one obtains another effective potential 

y m . ( H )m = f Y> + ^ (vj ^) ' + f + ^. (38) 

a \ a / 

The renormalized mass to order of (Ab^) 2 is given by 

2 _ d2 V(0) 



TO R(0) 



c=0 

2 A B ft 7 (A B fi) 2 7 7 
m B + ~rr J o 5 — JqJ-u (39) 



d 3 k 



where 

c 2 =m 2 +k 2 . (40) 
To compare Eq. @ with the GEP results, (Q and @, we first expand Eq. @ as 

0,(Mo) = «« + M _ J j_ x - %^J2 + 0(A B ^)3, (41) 



and then obtain 



4ui a 8ui a 32cj2 



^o(Mo) = Jo - Jo J-i + G(A B ?i) 2 , 



/-i(Mo) = J -i - + 0(X B h) 2 . (42) 



Hence, the renormalized mass 



ml = m| + ^ Jo - Jo J-i + <D(\ B h)\ (43) 

coincides with Eq. ( |39| ) up to (A?i) 2 . All higher order terms in the expansion of ml and Ar are a consequence of 
summing over daisy and superdaisy diagrams in GEP. 

Finally, we compare the LvN method with the time-dependent Hartree-Fock method. We divide the field and the 
momentum into a classical background field and a quantum fluctuation: 

$(t,x) = ^ c + $,(t,x), n(t,x) = 7r c + n 9 (t,x). (44) 

The Gaussian state has not only ($) = (f> c but also (II) = 7r c . The quantum fluctuations have the zero vacuum 
expectation value, = and = 0. In the Hartree-Fock method, $ 3 is replaced by 3($ 2 )<I > g . Then the field 
equation for quantum fluctuation is given by 

d»d>* 9 + (ml + + ^-m^j = 0, (45) 

and for classical field by 

d.d^c + Ul + ^(l» 2 ) J <f> c + ^</>l = 0. (46) 
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To remove the infinite term 

(*2> = ft5>*¥>„ = Wo(w>), (47) 

a. 

we add counter terms (8m 2 /2)<E> 2 and (5A/4!)<I> 4 , which is equivalent to writing the bare coupling parameters as 

m%=m^ + Sm 2 , Ab = Ar + £A. (48) 

By choosing the counter mass term as dm 2 = — (Ab7i/2)Io(mo)j we obtain from Eq. ( |45| ) the renormalized field 
equation for <p c = 

+ m£)$R, = 0. (49) 

Note that the renormalized mass obtained by adding the mass counter term is the same as Eq. ( p3|) in the GEP or the 
LvN method and the Fourier transform of Eq. (^) would be the operator counterpart of the renormalized mean-field 
equation 



III. ANALYSIS OF THE COUPLED QUARTIC OSCILLATOR MODEL USING MSPT METHOD 

The multiple scale perturbation theory (MSPT) was introduced to find solutions to nonlinear systems with a 
perturbatively small coupling constant, say Xh |25|,^6]]. The nonlinear systems exhibit distinct characteristic behaviour 
on different time scales determined by Xh. A conventional perturbation theory leads to secular (boundlessly growing) 
terms due to the resonant coupling between a lower order and the next leading order. However, the MSPT method 
provides a systematic technique of eliminating the secular terms and, as a consequence, gives rise to the solution 
with an effective frequency, which can be obtained by summing the most secular terms to all orders of Xh in the 
conventional perturbation theory J2f|. Even the lowest order MSPT solution is valid over a time scale l/(Xh). In this 
section we apply the MSPT method to analyze the system of two coupled quartic oscillators. It is a useful toy model 
for getting some insight into the phenomenon of mode mixing that one encounters in the study of nonequilibrium 
quantum fields during a spontaneous symmetry breaking phase transition. The MSPT was applied to a single classical 
as well as a quantum anharmonic oscillator in Rcf . p6| . In the latter case, the technique involves solving a complicated 
set of coupled operator equations. However, in the LvN method, the auxiliary equations for the Gaussian state or its 
excited states are the classical mean-field type equations or the linearized Hartree-Fock equations, which can be more 
easily tackled. As mentioned earlier (see also |2lj), these coupled equations are a direct consequence of the fact that 
the creation and annihilation operators given in Eq. (II) satisfy the LvN equation. 

The model to be studied in this section is described by the Hamiltonian 



H 




a 2 2 




(50) 



The Hamiltonian ( [50|) is obtained by taking two typical modes of scalar field from the field Hamiltonian The 
most physically interesting case is the coupling between the hard and soft modes. The approximate Gaussian state 
for the Hamiltonian (50) with (qi) = = ($2) is given by the product 



*o(<Zi,92) = #1 (91)^2(92) 



(51) 



of the Gaussian state for each mode: 



1/2 



{2<Khv* i v i yl i \v* 
The Vi, (i = 1, 2), satisfy the coupled equations 

Xh 



exp 



1 v 



i „2 



(52) 



Ul + WjVl + — (v\vi + V 2 V 2 )vi = 0, 
v 2 + ujIv 2 + -rr{v*Vi + V2V 2 )v 2 = 0, 



(53) 
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together with the Wronskian conditions 



Vi Vi - ViV: 



(54) 



Note that (Xh) in Eq. (|5^) is the nonlinear coupling constant. As each time scale of order {Xh) n t ~ 0(1) with 
n = 1, 2, • • • , shows a different characteristic behavior of the solution, we introduce a multiple of different time scales 



r (1) = (A?i)i, • ■ • , r (n) = (Aft) n f,- 



The solution to Eq. (53) is expanded in a series of Xh as 

Vi = (Xh^V^ (t, T(!) , ■ • • , T(„) ,•■■), 



(55) 
(56) 



n=0 



where is of the order of unity. The time derivative is now given by 
d 2 Vi _ B 2 vf ] t 2 B 2 vf ] ( 3t (1) \ t B 2 vf ] ( 9r (1) \, 2 



dV w (dr. 



(2) 



a< 2 a< 2 



dtdrm \ Bt 



Bt 2 
OT {i) 



Bt 



BtB 



T(2) 



t9/; 



, 2 3 2 ^ (1) f 9r (1) ^ { B 2 v^ ( 9r (1) ' 



at 2 0f07- ( i) i a< 



or (l) 



Bt 



9M 2) , 2 92 ^ (2) ( 9r (i) ^ , 9 \ (2 ] ( dT W ' 



at 2 



dtdr^i) \ Bt 



Bt 2 
OT (i) 



dt 



To zeroth order (Xh) , v { satisfies the simple oscillator equation 



•■(0) , 2 (0) n 



(57) 



(58) 



First, we consider the case wf > 0, which is the quantum mechanical analog of the coupling between two hard 
modes. As r = Tm is a time scale independent of t, one looks for the zeroth order solution to Eq. ( ^8| ) of the form 



=Ai(r)e iUit + Bi(T)e- iUit . 



Collecting terms of order (Xh) from Eq. (57), one obtains the nonlinear equation 



dV (0) 1 



.,(°K,(°) -L J°)*J°) \J°) { - 



- 9ld - + 2\ vi Vi +V2 V2 r 



= 0. 



(59) 



(60) 



The secular terms, which are proportional to e ; * in the second curly bracket of Eq. (|60f), should be eliminated to 
prevent v\ from growing boundlessly in time. This requires the coefficients to satisfy 





i 


Bt 




3B l 


—i 


Bt 


~ 4^^ 



B*B l+ J2 A*A 1+ B*B, U 

3=1,2 ) 



(61) 



3=1,2 



Equation ([H]) together with its complex conjugate for the coefficients A*(rm) and £?*(Tm), determines the initial 
values 



At^Mr) = A*(0)M0), Bt(r)Bi{T) = £,* (0)8,(0). 
Now the solutions to Eq. (pit) are given by 



(62) 



Ai(r) = Ai(0)exp 
Bi(r) = B l (0)exp 



— B* (0)8,(0) + ( A j(°) A M + B*(0)Bj(0)) 

1 I 3=1,2 



— ^(0)^(0) + £ W)^-(O) + S;(0)Bi(0)) 

1 I 3=1,2 



(63) 
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At order (A7i)°, the initial Gaussian state is given by 

Ai(0) = 0, Bi(0) 



Therefore, the solution for 



,(°) 



becomes 



(64) 



,(0) 



where each mode has a new shifted frequency 



a" = 



XTi x ^ 1 

Awi 2u>i 



(65) 



(66) 



J = l,2 



The second term can be obtained by summing the most secular terms to all orders in the conventional perturbation 
theory As each vf 3 ' satisfies the linear equation (58), we choose the coefficient to satisfy the Wronskian condition 
@ 

(67) 



(-) 



2ft, 

Note that the solution j67|), when expanded in power of Xh, yields the first order solution , as it has the correct 
frequency and coefficient up to that order. By comparing Eq. (p6J) with the gap equation 



ft,; = 



Xh 
T 



y — 

2ft 7 - 

j=l,2 3 



1/2 



(68) 



obtained from the exact solution to Eq. (| 



-iQ.it 



(69) 



it can be shown that oj - -* is the binomial expansion of ftj to order of A7L Hence the lowest order MSPT solution 
gives the correct frequency or the energy to order [Xh). 

We now consider the more interesting case of one hard and one soft mode. To mimic the interaction between the 
hard and soft modes after a sudden quench, one frequency squared, say cj|, will instantaneously change the sign from 
positive to negative at t — 0, i.e., u>\ — ► uj\ and wf — > 



§. Then Eq. (M) for t > takes the form 



2-2 2 

H = — + —qi + — 
2 2 yi 2 



2 9 2 2 + TT(9i+6g^ 2 2 + ^), 



4! 



(70) 



exactly the Hamiltonian for the interaction between a hard mode and a soft mode in $ 4 -theory which will be discussed 
in the next section. The equations for the auxiliary variables v\ and i>2 to 0((Xh) ) become 



"(0) . ~2~(0) 
V 1 + WiV{ 



0, 



~(°) ~2~(0) n 
v 2 - ^2 V 2 = °- 



(71) 



As for the former case of two hard modes, t = {XTi)t introduces another time scale in addition to t, and the solutions 
to Eq. (fnl) take the form 



i> { " ] = A(r)e iSlf + B{T)e' lQllt , 
4 0) =C , (r)e li2t + ^(r)e- i2t . 



(72) 



Secular terms appear in equations for due to the resonant coupling with vf^ . These are the terms proportional to 



e ±iu±t £ or ~( anc j e ±uj 2 t f or ^( anc j j eac j boundlessly growing solutions. By eliminating these terms, one obtains 
the equations 
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dA 


i 


97 




dB 


—i 








— i 






db 


1 


1h 


4tl> 2 



A* A + 2B*B + C*D + CD* jA, 
2A*A + B*B + C*D + CD* jf?, 
A* A + B*B + 2C*D + CD* jc, 
A* A + B*B + C*D + 2CD* }D. 



(73) 



A similar set of equations involving coefficients A*, B*, C* , D* , is obtained by requiring that the secular terms for 



v{ ' and v 2 



vanish, which is the complex conjugate of Eq. ([73|). Then it follows that 

A*{t)A{t) = l*(0)i(0), B*(t)B(t) = B*(0)fl(0), 
C*(t)D(t) = C* (0)5(0), <5(t).D*(t) = C*(0)L>*(0). 



(74) 



As the system evolves from an initial Gaussian state for t < to another Gaussian state for t > after the 
phase transition, the Gaussian state (|52|) should be continuous across t = 0. The continuity of the Gaussian state is 
equivalent to the continuity of Vi and Uj at i = 0, which determines the coefficients of the zeroth order solutions after 
the phase transition in terms of the initial data: 




2^2U} 2 \ Ll>2 

By substituting Eq. (|75| ) into Eq. ( |T4D and solving Eq 

»-;( )/ 



«w(t) = A(0)e lU i * + B(0)e 
^\t) = C(0)e^ h + D(0)e- il 



2y/2uj 2 \ W2 

, one finally obtains 



(75) 



Q\~'t 



(76) 



where 



(±) _ ,-, 




n 



By requiring the Wronskian conditions ( pif ) and matching the first order solutions (| 
coefficients as 



(77) 

we find the correct first order 



Si(t) = 



1 



29. 



(-) \ n 



29 



(-) 











9^- 






-n<->. 



1 



2ft 



(-) 



2ft 



(-) 



fti + ^ 




ft«H 


-**$-> J 


9^- 


-ift[ _) ' 







(78) 



A few comments on the interpretation of solutions (|65|) and ( |78j ) are in order. These solutions are valid till 
t ~ l/(Xh). The coupling between the hard modes just modifies the frequencies ( |66[ ) even after the phase transition. 
However, the coupling between the hard and soft modes after the phase transition exhibits some interesting features 
for the both modes. One interesting feature is the disparity between the positive and negative frequency given by 
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5£li = Q 



(+) 



f2 



A/8tD 2 for the hard mode and SQ2 = ^4 + ' > ~^ 



>(-) 



A/8cD| for the soft mode. 



This difference is 



a consequence of the resonant coupling for the positive and the negative frequency solution and for the exponentially 
growing and decaying solutions. Another interesting feature is the appearance of an imaginary part of Q 2 f° r the 
soft mode. This implies that V2 begins to oscillate due to the coupling between the hard and hard modes. The time 
scale of oscillation is estimated to be (16w|)/(A?i), which is larger by the factor 16w| than the time scale 1/Xh for the 
lowest order solution to be valid. The period of the oscillation is too long to be observed. 

Though the time scale for oscillation is outside the validity region of the zeroth order solution, the oscillatory 
behavior gives more information about the later stage of evolution towards and oscillation around the true vacuum. 
The numerical analysis of the first order solution shows that the early stage of exponential growth due to instability 
is eventually terminated by the back-reaction from nonlinear terms pi] ]. It is the oscillatory part that determines 
the absolute width of the Gaussian state after the field samples the true vacuum. To see how the oscillatory part to 
restrict the Gaussian state, we use the lowest order solution u 2 , which is the auxiliary equation for the Gaussian state 
for the soft mode. The correct physical interpretation may follow from the Gaussian state (p2n: 



1 



(2tt^*C 2 ) 1 / 4 I v. 



- \ 1/2 

ex p 



A 



32<D| 



2^2 

T10J2 



^Re(o4 +) )e- 2Rc ^ +) )* 



+j{5 ^ (1 



3 -4Rc(ii^ +, )t 



(79) 



Immediately after the phase transition, the width (dispersion) of the Gaussian state (|7^) is largely determined by the 
second term of the real part which originates from the exponentially growing part of V2- But, as the phase transition 
continues, the first term from the imaginary part of f2 2 i s comparable to the second term when 



-2Rc(<i< +) )t 



Xh 

64wf 



(80) 



The time scale of Eq. ( |80| ) is outside the validity of the lowest order solution. Before reaching the time 
the back-reaction (Xh/2)(v2V2)v2 of the self-interaction begins to dominate over — a)|«2- So the more accurate first 
order solution oscillates around the true vacuum after an exponentially growing period and the Gaussian state has 
approximately a constant width, like the first term of the real part. That is, the Gaussian state stops spreading. 



IV. MSPT AND RENORMALIZED FREQUENCY 



The main advantages of the MSPT method, especially within the framework of the LvN formalism, was discussed 
in the previous section. Nevertheless, it is necessary to make some comments on this method in the field theoretic 
context. The process of eliminating the secular terms at various orders of the MSPT method is equivalent to summing 
certain diagrams to all orders in the conventional perturbation theory. It is precisely because of this fact that useful 
information pertaining to the divergence structure of the theory can be extracted even from the lowest order solution 
in the MSPT method. As will be shown below, the lowest order MSPT solution of a quantum field yields a shifted 
(effective) frequency which contains divergent contributions and the removal of the most divergent term essentially 
amounts to mass renormalization. 

The requirement that the creation and annihilation operators (11) satisfy the LvN equation, leads to the mean- field 
equation of the auxiliary field variable of each mode 




where 

^a=™l + k 2 . (82) 

The exact solution to Eq. (^l|) for the static system without phase transition is given by Eq. ( |l6| ) . Now we apply the 
MSPT to solve perturbatively Eq. (^l]) in a power series of Ab?i and to find the correct renormalized frequency. As 
in Sec. Ill, introduce a multiple of time scales 

r (1) = {X B h)t, ■ ■ • , r (n) = (X B h) n t, • • • , (83) 
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and expand the solution in the power series of Xh 

fa = $3(ABa)VS' ) (*.T(l)'" , ' 7 '(»)' ,, 0- 



As we shall find the solution (|8J) only to the first order, we denote T(i) = r and truncate it to that order 



f a (t) = f^(t) + (X B fi)f^(t,r). 



The time derivative is now given by 



dt 2 dt 2 dtdr [ dt 



(A B ») 



a 2 ^ 
dt 2 



The terms to order (AbS)° satisfy the equation for a simple oscillator 

The general solution is a superposition of the positive and negative frequency solutions 

<pW =A a (r)e i ^ t + B a (T)e- i ^ t . 
At the next order the terms of order (Ab?i) in Eq. ( [sT] ) obey the equation 



° fa 2 (1) 



d 2 ^ 1 

<9i<9r 2 



+ ^(E4 0) v?>i 0) 



0. 



(84) 

(85) 
(86) 

(87) 
(88) 

(89) 



As ifa ^ has the same natural frequency u> a as fa\ any term proportional to e ±l "°* in the second curly bracket is a 
secular term to fa . To avoid unnecessary growing terms, these secular terms should be eliminated from the second 
curly bracket. This requires 



dA a 


i 


dr 




dB a 


—i 


dr 


4cj a 



B a B a + J2A*pAp + B*pBp \A a , 



- = j-{ A* a A a + A* Ap + B* p B p \B a . 



As the coefficients of A a and B a are real, Eq. ( pfj| ) is unitary, so has the solution 

A* a (r)A a (r) = A* (0)A a (0), B* a (r)B a (r) = B*(0)fl a (0). 
Hence the solutions are given by 



(90) 



(91) 



A a (t) = A Q (0)exp 
B a {r) = B a (0) exp 



_|S*(0)B Q (0) + ^^(0)^(0) + ^(0)^(0) 



4w Q 



A* (0)A Q (0) + J] ^(°)^(0) + B* fi (Q)B p (Q) 



As the Gaussian state before the phase transition has the initial data 

1 



Aa(0) = 0, B a (0) = 



the solution becomes 



(92) 



(93) 



(94) 
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where 



i / , w Q + - — Jo- 



By requiring the solution fl94J) to satisfy the Wronskian condition d27|) , the solution is given by 

1 



2ft 



(-) 



: exp[— ifi^ 'i] 



(95) 



(96) 



The solution ( p^ ) of the MSPT method coincides with the frequency ( p4"| ) to the first order of XrTi. As a consequence, 
the MSPT solution (|9^) recovers the GEP ( |l8| ) and renormalized mass ( |23| ) to the first order Ar?i. All higher order 
corrections to the MSPT solution are expected to exactly recover the solution ([lfj]) of GEP. In this sense the MSPT 
method can be understood as a powerful tool to systematically solve Eq. (|2^ ) even for the nonequilibrium case with 
time-dependent coupling parameters. 

We now apply the MSPT to a quantum field changing suddenly from one parameter m B to another m B > at 
t = 0. The mean field equation (|8l]) still applies with the modification 



and the solution 



has the form 



~ 2 ~ 2 i i 2 



^°)(t) = l ce (r)e^ t + J B Q (r)e- 



The solution is given by 

A a {r) = i Q (0)exp 

B a {r) =B Q (0)exp 



(0)fl a (0) + A* (0)A fi (0) + B% (O)B (0) 



g | i* (0)i Q (0) + £ A% (0)A P (0) + (0)^ (0) 



The initial data is determined by the continuity of tp a and (p a at t = 0: 



A a (0) 



1 



1-^ , Sa(0) 



1 



2V2t^ 



1 



The zeroth order solution can then written be in terms of the modified frequencies as 

^\t) = A a (0) exp {in { +h) + B a (Q) exp[-iQ^h]. 

where 



Q(±) _ ^ , M 



\^ 1 (to| - to|) 

1.0 2 



The shifted frequency (|102|) can be written as 



q(±) _ s , M 



where we used the relation 



(m| - ro|) (^r~to|) : 
J o H ^ J-i H ^ 



.y 1 +J_fi ± ^Y 



™2 



(97) 



(99) 



(100) 



(101) 



(102) 



m B — m B = ?ti r — rn 



(103) 



(104) 



and rewrote the sum 
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1 f$ - Wjl 



/3 p P /3 P /3 



J_ 1 + (in >-^)£ * . 

/3 P P 



(105) 



Finally, by fixing the coefficients of Eq. ( 101 ) to match Eq. (|96| ) and satisfy Eq. (|27| ) we obtain the lowest order 
solutions 



1 



20 



(+) , o(-) 



(-) vn^ ; + fi£ 



exp[ifi(+>t] 



1 



(+) , o(-) 



(-) in^ + ni," 



exp 



(106) 



From the GEP calculation, Eqs. ( |19| ) and (|35|), of Sec. II , it follows that the new renormalized frequencies can be 
written as 



(ro| - m|) (m|-m|) 2 

2 _1 2 ^ ' 8w Q 



x - 1 _L_ 

2^ 0,,,3, r ,2 + Q. , 



1± 



(107) 



V. RENORMALIZED SOLUTION AFTER PHASE TRANSITION 



After the quenched phase transition at t = 0, the soft modes (long wavelength) with momenta k < % become 
unstable and evolve out of equilibrium, whereas the hard modes (short wavelength) with k > evolve towards a 
new equilibrium. The Gaussian state after phase transition can then be determined by two sets of the auxiliary fields, 
one for the stable hard modes and the other for the unstable soft modes. The mean field equation ( |8l| ) for the hard 
modes then become 




where 



and for the soft modes 



where 



^Ha + u Ha m« + — [ 2^ vptppj m a = o, (108) 



^ a =^-mi, (109) 



VSa-Us a VS a + ^(^&0&l3j&S a = O, (HO) 



u>l a = m| - k 2 . (Ill) 



Once again, by introducing two independent time scales t and r in accordance with the MSPT method, the solutions 
of zero order 0((Xb^) ) for the hard and soft modes can be written as 

$1 = Mr)e 1 ^ 1 + B a {T)e-^«\ (112) 
4°2 = ^(r)^" 1 + D^e-*^. (113) 



Substituting these zeroth order solutions into Eqs. (108) and (|ll0] ) and eliminating the secular terms e ±MH °* for (p^ a 
and e ±WScl * for <fa^, one finds the equations for coefficients: 
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dA a 



OT 


4WHa 


dB a 


—i 


OT 




dC a 




dr 


4wSa 


db a 




dr 


4wSa 



KB a + 53 i;^ + + ^ qs^ + Cpfy \A a , 

0> 0< ) 

A* a A a + 53 A^Ap + b;b p + 53 C* p Dp + Cpb* U a , 



- i c;£ Q + 53 A}Ap + b;b p + 53 c;b + c f3 b; \c a , 



— \ D* a c a + 53 A} 3 Ap + Bpbfj + 53 c;b p + c b% \b 



(114) 



0< 



Here /?> and /3< denote the restricted sum (integral) over the hard and soft modes, respectively. Another set of 
equations, the complex conjugate of Eq. ( ]114| ), can be obtained for the coefficients A* a , £?*, (7* and D* by eliminating 
the secular terms for <p^/ and V>g„ . This set of equations together with Eq. (114) determines the initial values of 
the coefficients: 



A* a (r)A a (r) = i* (0)A a (0), B* a (r)B a (r) = B* a (0)B a (0), 
C* a (r)b a (r) = C£(0)A,(0), C a (r)D*(r) = C a (0)D*(fl). 



(115) 



A a (r) = i Q (0)exp 
B a (r) = B Q (0)exp 
C„(r) = C a (0)exp 
£>a(r) = A*(0)exp 



By using the above relations, we find the solution of Eq. (114) as 

B* a B a + 53 A}A P + B* B + 53 C* p b p + C D* 1 

0> 0< ) 

A* a A a + 53 A%A + B* p B + 53 C* p Dp + CpD* \ 
0> 0< ) 

- 1 c*b a + 53 A;A f3 + b;b + 53 c* b + c p b% 1 

- J ^;c Q + 53 A* A P + b;b p + 53 c;b p + c p b* 1 
Sa I 0> 0< J 



4^H Q 



4^H Q 



4ws c 



(116) 



The initial data (115) are determined by continuity of (fia^ and 

1 



A a (0) = 
C tt (0) = 



1 



2v2cj 

, \ ,-, 1 / , . . LOa 

— ; 1 + l ~ — 

2V2w Q \ uj Sc 



, A,(0) = 



(117) 



It is worth noting that the integral for the soft modes, the last two terms of the summation in the exponent of Eq. 
(116), is a restricted one with its upper limit being tob- On the other hand, the integral for the soft modes, the first 
two terms of the summation in the exponent of Eq. (116), has me as its lower limit of integration. For k 3> %; onc 
can write wh/3 ~ u>/3, whence the ultraviolet (quadratic) divergence mainly comes from f° r the hard modes. 

In terms of new frequencies for the hard modes and for the soft modes, the zeroth order solution can be 
written as 



0Ha = ^(0)e in H.'* + B a (0)e- m *° 



4°2 =C Q (0)e s ^' t + D Q (0)e^ 



>t 



(118) 



where 
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o(±) - r, 



r>(±) 



Ei + 



(m| + m|) 



4(2' 



So 



2 W/3 2 



\ - 1 (m|+m|) 
^ 2ujr 2 



E 



1 



2u; 



y 1 

3 2w /3^S/3 



E 



2w 



1 

1 

8oj q 



1 ± 



2-i 



2n 



1± j- 



(119) 



Here ^2 a denotes the restricted integral over soft modes, having a finite upper momentum cut-off given by ro B and, 
therefore, making a finite contribution to the renormalized frequencies. However, the quadratic ultraviolet divergent 
term has to be appropriately regularized to yield the renormalized frequencies. Contrary to a naive belief that the 
frequency of the soft modes may be finite, the soft frequency also has the same ultraviolet divergence as the hard 
modes and renormalization is required to remove it. This is a consequence of the nonlinear coupling between the 



soft and the hard modes. As the mass counter term dm 2 has the same sign, 



5m , even after the phase 



transition, the sum of bare squared masses after the phase transition equals to that of renormalized squared masses: 

(120) 



2 i ~ 2 



2,-2 
m R + m R- 



Using the identity 



E 
E 



2u)pCj 



~,2 
H/3 



2u!ffLU. 



~2 
S/3 



0> 2 ^ 

-y— 



ml 



i 



2lo\Cj? 



(121) 



we write the shifted frequency as 



n(±) - ,~, 4- Ab ^ 



5 2u}3 ^ 



n(±) 



(m| + m|) 
•Jo H ^ J-i 



(m£ + ) 2 



E 



i 



/3> 2 <^H/3 



E 



i 



1 

8WQ, 
8w Q \ 



1 ± 



l±i- 



^Sa 



(122) 



The first two terms of the MSPT frequency ( |122| ) 



4w„ 



(123) 



gives, to the first order of XbTi, the renormalized frequency 



h 2 



A B 7 i 
~2~ 



-2 , '- T 

U~ ± Jo- 



(124) 



Finally, the lowest order solutions which match the solutions (|96|) before phase transition and satisfy Eq. 
given by 



arc 



<PHc 



1 



Ha Q 



2n ( a ) l^Ha 



(-) 



20 Q } V ^ 



(-) 

Sa 



1 



n 



(+) 

Ha 







(")' 



2fii } V^Ito 



(-) 



n 



"Sa 



n 



(-) 

Sa 



(125) 



It is shown that the MSPT method, when applied to the mean field equation (81), yields the correct renormalized 
frequency even after the phase transition. 
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VI. CONCLUSION 



In summary, we have studied certain aspects of the renormalization of quantum helds evolving out of equilibrium 
using the recently developed Liouville-von Neumann method [ pi) . The multiple scale perturbation theory (MSPT) 
provides a powerful and useful tool in finding an analytic approximate solution to nonlinear field equations. The mean- 
field equations for the auxiliary variables for the Gaussian state, which are equivalent to the vacuum expectation of the 
Hartree-Fock equation, are analyzed using the MSPT. As a simple quantum mechanical model, two coupled quartic 
oscillators are studied. In particular, the model with one unstable quartic oscillator, mimicking a soft mode of a 
nonequilibrium field, sheds some light on the effect of nonlinear mode coupling between the soft and the hard modes 
in the nonequilibrium evolution of phase transitions jnj. The MSPT method is employed to the $ 4 -theory before the 
phase transition and its lowest order solutions to the mean-field equations involve the correct renormalized frequency 
consistent with the Gaussian effective potential |^3[ . Also the MSPT method turns out to be a powerful and practical 
tool for the renormalization of nonequilibrium quantum fields such as they appear in phase transitions. 

Although we have used Eq. (|2^) in order to express the modified frequencies in terms of the finite renormalized 
parameters, there is a caveat which is necessary to keep in mind. The renormalized mass is obtained from the 
quantum fluctuations of the Gaussian state around 4> c — for the static and stable <I> 4 -theory. When the symmetry is 
spontaneously broken, (f> c — is no longer the true vacuum state. However, the justification of Eq. (|23|) in obtaining 
renormalized parameters, is solely dependent on the same Gaussian approximation whose quantum state is the time- 
dependent Gaussian state around <p c = throughout the phase transition. This effectively amounts to studying the 
nonequilibrium field dynamics during the rolling down of the field and before the fluctuations start probing the true 
vacuum state. As far as the soft modes are concerned, they just see an inverted potential. Therefore, the lowest order 
solution of the MSPT method exponentially grows indefinitely and begins to oscillate slowly. The imaginary part of 
the solution, leading to the oscillation, determines the width of the Gaussian state at later times, roughly covering 
the true vacuum state. The exponential growth of the lowest order solution is terminated by the first order correction 
of the MSPT method when the field begins to sample around the true vacuum state |n| . 

In this paper we confined our attention to the Gaussian approximation using the LvN formalism (which is equiv- 
alent to the Hartree-Fock or mean-field method) to address issues related to renormalization of a spontaneously 
broken $ 4 -theory. We now point out some of the limitations of the Hartree-Fock method since an understanding of 
those limitations is crucial in developing techniques for more accurately describing nonequilibrium dynamics of phase 
transitions. As has been recognized in the literature, the Hartree-Fock method (sometimes called the collisionless 
Hartree-Fock method) neglects the effects of scattering. It also cannot account for late time thermalization of the 
system since it takes into account the effect of interactions only through a mean field. It is based on the assumption 
that scattering effects do not play a very crucial role in the early stages of phase transition dynamics. However, recent 
studies [^H indicate that scattering effects may be crucial to understanding the dynamical evolution of the system far 
from equilibrium. Similarly, in the Liouville-von Neumann method, the quartic part Hp, which is a perturbation to 
the quadratic (Gaussian) part Hq and is negligible before the phase transition, grows comparable to Hq and signifi- 
cantly contributes to the Gaussian state |33|j . These terms may correspond to the next-to-leading order contribution 
to the effective action of Ref. p2[ . The study of non-equilibrium phase transition dynamics beyond the Gaussian 
approximation is under investigation. 
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